A one dimensional hard-point gas as a thermoelectric engine 
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We demonstrate the possibility to build a thermoelectric engine using a one dimensional gas of 
molecules with unequal masses and hard-point interaction. Most importantly, we show that the 
efficiency of this engine is determined by a new parameter YT which is different from the well 
known figure of merit ZT. Even though the efficiency of this particular model is low, our results 
shed new light on the problem and open the possibility to build efficient thermoelectric engines. 

PACS numbers: 74.25.Fy,84.60.Rb,05.60.Cd,05.45.Pq 



A combination of mathematical results, numerical 
studies and even laboratory experiments have greatly im- 
proved our understanding of phenomenological transport 
equations in recent years. The derivation of such equa- 
tions from purely dynamical laws, classical or quantum, 
has been one of the main subject of interest [l|]. Even 
though a complete rigorous picture is still lacking, it is 
however clear that dynamical chaos is an essential ingre- 
dient. 

A better understanding of the above problem is impor- 
tant not only from a fundamental point of view in order 
to provide a justification of phenomenological laws. It 
is also relevant for several applications. One particular 
important aspect is the connection with thermoelectric 
power generation and refrigeration @, Indeed, due 

to the increasing environmental concern and energy de- 
mand, thermoelectric phenomena are expected to play an 
increasingly important role in meeting the energy chal- 
lenge of the future. However, the main difficulty is the 
poor efficiency of existing devices. Indeed the suitability 
of a thermoelectric material for energy conversion or elec- 
tronic refrigeration is characterized by the thermoelectric 
figure of merit Z = aS 2 /k, where a is the coefficient of 
electric conductivity, S is the Seebeck coefficient and k is 
the thermal conductivity. The Seebeck coefficient S, also 
called thermopower, is a measure of the magnitude of an 
induced thermoelectric voltage in response to a temper- 
ature difference. 

For a given material, and a pair of temperatures Th 
and Tq of hot and cold thermal baths respectively, Z is 
related to the efficiency rj of converting the heat current 
Jq (between the baths) into the electric power P which 
is generated by attaching a thermoelectric element to an 
optimal Ohmic impedance. Namely, in the linear regime: 
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where ?7carnot = 1 — Tc/Th is the Carnot efficiency and 
T = (Th + Tc)/2. Thus a good thermoelectric device 
is characterized by a large value of the non-dimensional 
figure of merit ZT. However, in spite of the fact that 
the second principle of thermodynamics does not impose 
any restriction on the value of ZT, all attempts to find 
high ZT values (let us say ZT > 3 at room temperature) 
have failed. The problem is that the different transport 
coefficients S, a and k are interdependent, making the 
optimization of ZT extremely difficult. We believe that 
a better understanding of the possible microscopic mech- 
anisms [5] which determine the value of ZT may lead to 
a substantial improvement. 

In a recent paper Q , a dynamical system approach to 
a Lorenz gas type model has been used and an interesting 
mechanism to reach high ZT value has been discovered. 
More recently, a one-dimensional, di-atomic disordered 
chain of hard-point elastic particles has been considered 
[7J and it has been found that ZT diverges to infinity with 
increasing the number of particles inside the chain. This 
result suggests the possibility to build a thermoelectric 
engine by connecting two heat baths with two chains of 
different sizes. 

In this paper we analyze such an engine and we show 
that, indeed, a non zero circular current is established 
inside the system when the stationary state is reached. 
Quite interestingly the efficiency of such an engine ap- 
pears to be unrelated to the figure of merit ZT and it 
remains quite low in spite of the fact that ZT becomes 
larger and larger with increasing the system size. Indeed 
we analytically show and numerically confirm that the 
efficiency of such engine is determined by a new figure of 
merit YT. From one hand the low value of YT explains 
the poor efficiency of the engine; on the other hand, it 
sheds new light on future directions for increasing the 
efficiency of thermoelectric engines. 

We consider a simple model of the engine (see Fig. 
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FIG. 1: Scheme of the thermoelectric engine studied in this 
paper. It consists of two heat baths at different temperature 
Tl and Tr and two channels of different lengths. See text for 
the details. 



[T]) which consists of two heat baths and two connecting 
channels A and B of length La and L b respectively. We 
assume that the whole system does not exchange parti- 
cles with the outside. Therefore the total number of par- 
ticles N in the two channels and in the two baths remains 
constant during the simulation. The two heat baths are 
kept at different temperatures Tl and Tr, and have the 
same finite length V. Each channel is a one-dimensional, 
di-atomic disordered chain, of hard-point elastic particles 
with velocities Vk and masses rrik € {Mi,M%} randomly 
distributed (we use convenient non-dimensional units). 
The particles interact among each other through elas- 
tic collisions only. A collision between two neighboring 
particles with mass rrik and mi causes a change of their 
velocities Vk and?;/ as Avk = 2mi(vi—Vk)/(nik + mi) and 
Avi = 2mk{vk — vi)/(m,k + mi). An efficient algorithm 
has been developed which performs correct chronologi- 
cal order of collisions and update particle's positions and 
velocities in ~ log N computer operations per collision. 
We fix the length of the channel B to be small and study 
how the properties of the engine vary as the length of 
the channel A is increased. In particular, if channel B is 
small enough then particles can pass trough it without 
suffering any collisions. 

The heat baths are modeled in the following way. Sup- 
pose n\ particles of mass M\ and ni particles of mass M2 
are confined in a box of length V. The particles collision 
rate with one end of the box is given by 
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where 7* = pJ Pi = n i/ V (* = !> 2 )> P = n / v and 
n = Ti\ + Hi- Then the heat bath model is straightfor- 
ward: if one end of the box is opened, the particles are 
emitted with the same rate 7, and the time interval t 
between two consecutive emissions is a random variable 
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with to — I/7. The mass of the emitted particle is as- 
signed to be Mi (i = 1,2) randomly according to the 
probability 
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and its velocity is generated from the distribution 
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As expected, when M\ — M2, the heat bath model of 
identical particles is recovered. 

The emission rate from the left heat bath into channel 
A is therefore jl = ^^7i,z, + ^7^72, L where pl is the 
total particle number density at the left heat bath, pi^ 

is that of particles with mass M$, and j^l = Pl\J j^itr 
(i = 1,2). Similarly one can write the expression for the 
emission rate 7_r from the right heat bath. 

As to the channel B we assume that the emission rates 
are proportional to 7L and jr: 



7f 



(6) 



where r is an adjustable parameter. For r — channel B 
is in fact closed hence there is no net particle current. 

Conversely, whenever a particle from each channel ar- 
rives at the border with the bath it is absorbed by the 
bath, so that a stationary state is established after suf- 
ficiently long time. If the stationary state is such that 
there is a net particle current around the system, then 
one can use it to extract work. To this end we insert 
inside the channel B an auxiliary potential Ub which the 
particles have to climb thus performing useful work. Let 
us first consider the case where the channel B is short 
enough (i.e. Lb = 1) so that the probability of particles 
collision inside it is negligible. Suppose Tl > Tr then the 
current runs clockwise at Ub = due to the pressure bal- 
ance between the two heat baths. Thus a particle emitted 
from the right heat bath can pass through the channel B 
only when its kinetic energy is larger than Ub, otherwise 
it will return back. Therefore the probability with which 
the emitted particle passes through the channel B is 
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for both particles of masses Mi and M%. Hence the work 
extracted per unit of time is 



P = t^rVrl - 1l)U B - 
The efficiency of the engine thus reads 
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FIG. 2: Efficiency of the engine model illustrated in Fig. ((TJ) 
as a function of the potential Ub for La=4, 8, 16, 32, 64, 100, 
200 and 400 (from below to up). T L = 1.2, T r = 0.8, and 
r — 1. The length of channel B is unity. 

where Jq , and Jq , are the net thermal energy flows from 
the left heat bath in a unit time into the channel A, and 
B, respectively, which can be measured by numerical sim- 
ulations. Notice that, in particular, Jq can be expressed 
analytically through 

Jq.b = r(T L -f L - VrlTrJr). (10) 

As it is seen from Fig. [2] our engine works. Here we 
plot the numerically computed efficiency rj against the 
potential Ub for various lengths of the channel A. The 
solid curves are polynomial (3-order) fittings of the nu- 
merical data. 

In this calculation the length of channel B is set to 
unity and the length of both heat baths is very large, 
namely V = 10 3 . The total number of particles in the 
system is set to be N = 2V + La + Lb such that the over- 
all averaged particle density is unity. The two types of 
molecules are set to be equal in number and have masses 
Mi = 1 and M2 ~ 0.618 respectively. In addition, for 
a given length La of the channel A, there is an optimal 
value of Ub such that the efficiency is maximized. De- 
noting the maximum efficiency by 7? max , it is interesting 
to compute the dependence of ?7 max on La and compare 
it with theoretical expectations. 

Extensive and accurate numerical computations, sum- 
marized in Fig. [31 show that even though ry max increases 
with La, i.e. with the average number of molecules in the 
channel A, the increasing rate slows down very fast. Also 
it seems hard to find a simple fitting for the dependence 
of r) max on L A - 

We should notice that r] max depends also on Tl and 
Tr. However if 77 max is rescaled to the Carnot efficiency, 
its behavior with La is the same for all (Tj,, Tr) pairs, 
given that the temperature difference is small enough. 
This can be seen clearly in Fig. [3) 

Remarks: We have also investigated carefully several 
variants of the engine model. First, we have studied the 
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FIG. 3: The maximum efficiency of the engine for different 
temperature gradients. Here r—1 and Lb = 1- The solid 
dots are obtained from Eq. (|14|l and Eq. I)15[l. where the 
parameters k, a, S for both channels are numerically obtained 
as in Fig. Q. As it is seen there is perfect agreement with 
the directly, numerically computed, efficiency. 

dependence of 77 max on r. We have found that indeed 
??m ax can be slightly improved by adjusting r (e.g. ~ 30% 
as compared to Fig. however the efficiency remains 
very low. Then we have checked whether adding a po- 
tential U a against the current in the channel A can im- 
prove efficiency. In this case a particle in the channel A 
would undergo a parabolic motion between two consecu- 
tive collisions with its neighbors. This variant does not 
increase the maximum efficiency ?7 max either which has 
been found to depend only on the sum U a + Ub rather 
than on Ua and Ub separately. Also a longer length of 
channel B has been considered but it turns out that for 
a given La + Lb, the efficiency reaches its highest value 
when Lb < 1. To summarize, the various modifications 
to the engine model we have considered do not improve 
the efficiency significantly. 

Let us now explain why our engine works at a rather 
low efficiency in spite of the fact that ZT goes to infin- 
ity. In order to derive a theoretical expression for the ef- 
ficiency we use thermodynamic linear response relations 
for the heat currents Jq.a, Jq,b, and for the particle cur- 
rents J Pj a, Jp.B, in the channels A and B, respectively, 
namely 

Jq,a = -k a AT/La-T<7 A Sa(Ah/La + Ua) 
J p ,a = -cfaSa&T/La - <ja(A[J>/La + Ua) 

Jq,b = -k'b&T/Lb-T<j b Sb{&h/Lb-Ub) (11) 
J P ,B = -(tbSb&T/Lb - <tb(A/i/L b — Ub) 

where k' a = ka + TuaS\, k' b — kb + TobS 2 b , [i is the 
chemical potential, and AT = Tr — Tj,, A/i = fiR — /i^- 
In addition, in the stationary state we must also have 

Jp.A + Jp.B = 0. (12) 

The relations given by Eq. (fTTj) and (fT2")l represent 
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FIG. 4: Dependence on the channel length of various param- 
eters of the ID mixed gas. From these results it is clear why 
when ZT can be very large for a long channel, the new merit 
YT that governs the efficiency of the engine model (Fig. [TJ 
takes a lower value instead. 



five linear equations in the five unknowns Jq,a, Jq,b, 
Jp,Ai Jp,B, and A/x, whereas the other parameters, such 
as transport or Onsager coefficients are considered to be 
known. From the solutions of the above equations we 
compute the efficiency 77 as a function of the two bias 
potentials: 



r)(U A ,U B ) = — 



J p ,aUa — J p ,bUb 
Jq,a + Jq.b 



(13) 



Since r](U A ,U B ) is a quadratic function, we can easily 
maximize it by computing the optimal values of U A , U b . 
The final result can be expressed in a simple form as 



?7max/?7Car„ot - 1 - 2y/ (YT)~* + (YT) -1 + 2(YT)- 1 

(14) 

with a new figure of merit 



YT = 
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From Eq. (fT4"|) it can be seen that efficiency is given by 
this new figure of merit exclusively and it has nothing to 
do with the merit ZT. Clearly, as YT — > (YT -> oo) 
we have ?7 max -> (?7 max -> rfcamot)- 

In order to understand why our engine model has a 
low efficiency, we investigate carefully the dependence of 
the parameters a, k and S on the channel length. These 



quantities can be obtained from the Onsager coefficients 
[6j which in turn are calculated by measuring the par- 
ticles current J p , and the energy current J u , when the 
stationary state of the system is reached through 



J u = L uu Af3/L- L up Aa/L 
J p = L pu Af3/L - L pp Aa/L 



(16) 



with a — p/T, (3 = 1/T and p the chemical potential: 
p/T = (pilnpi + p 2 lnp 2 )/p - InVT + const. (17) 

In Fig. [4] we present the numerical results for a ID 
mixed gas chain with M\ = l,M 2 « 0.618, T = 1 and 
2pi = 2p2 = p = 1. The constant in Eq. (fT7|) is set 
to be ln2 + 1 such that p — 1 . The results clearly show 
that, as the channel length is increased, while the ratio 
a I k increases as a power law, the Seebeck coefficient S 
undergoes a slower and slower increase. As a result, while 
ZT increases indefinitely with the channel length L, the 
merit YT remains small. 

In summary, for the engine model (Fig. [T]) based on a 
one dimensional mixed gas, the efficiency is governed by 
a new figure of merit. Based on its relation to the ther- 
moelectric parameters (Eq. (fl"5j)). it can be expected that 
an efficient mixed gas type engine should be such that its 
Seebeck coefficient changes fast when the channel length 
is changed. Such a mixed gas might be that consisting of 
three, or more, types of molecules. We believe that our 
method of analysis of thermoelectric or thcrmochemical 
heat engines should be applicable to a wide range of mod- 
els which consist of two transport channels between a pair 
of baths. 
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